Tail-induced spin-orbit effect in the gravitational radiation of compact binaries 
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Gravitational waves contain tail effects which are due to the back-scattering of linear waves in 
the curved space-time geometry around the source. In this paper we improve the knowledge and 
accuracy of the two-body inspiraling post-Newtonian (PN) dynamics and gravitational-wave signal 
by computing the spin-orbit terms induced by tail effects. Notably, we derive those terms at 3PN 
order in the gravitational-wave energy flux, and 2.5PN and 3PN orders in the wave polarizations. 
This is then used to derive the spin-orbit tail effects in the phasing through 3PN order. Our results 
can be employed to carry out more accurate comparisons with numerical-relativity simulations and 
to improve the accuracy of analytical templates aimed at describing the whole process of inspiral, 
merger and ringdown. 
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I. INTRODUCTION 



A. Motivation 



During the last ten years a network of ground- 
based laser-interferometer gravitational-wave detectors 
has been built and has taken data at design sensitivity. 
It is a worldwide network composed of the Laser Interfer- 
ometer Gravitational wave Observatory (LIGO), Virgo, 
GEO-600, and TAMA and it has operated in the fre- 
quency range 10-10'^ Hz. Coalescing binary systems com- 
posed of black holes and/or neutron stars are among the 
most promising sources for those detectors. By 2016 the 
gravitational- wave detectors will be upgraded to a sensi- 
tivity such that event rates for coalescing binary systems 
will increase by approximately a factor one thousand, 
thus making likely the first detection of gravitational 
waves from those systems. In the future, space-based de- 
tectors like LISA should detect supermassive black-hole 
binary systems in the low frequency band 10~^-10^^ Hz. 

The search for gravitational waves from coalescing bi- 
nary systems and the extraction of source parameters 
are based on the matched-filtering technique, which re- 
quires a rather accurate knowledge of the waveform of 
the incoming signal. In particular, the detection and 
subsequent data analysis are made by using a bank of 
templates modeling the gravitational wave emitted by 
the source. The need of a faithful template bank has 
driven the development of accurate templates over the 
last thirty years. 

The post-Newtonian (PN) expansion is the most 
powerful approximation scheme in analytical relativ- 
ity capable of describing the two-body dynamics and 
gravitational-wave emission of inspiraling compact bi- 
nary systems [1]. The PN expansion is an expansion 
in the ratio of the characteristic orbital velocity of the 
binary v to the speed of light c. However, as the black 
holes approach each other toward merger, we expect the 



PN expansion to lose accuracy because the velocity of the 
holes approaches the speed of light. At that point, nu- 
merical relativity (NR) [2-4] plays a crucial role providing 
us with the dynamics and gravitational-wave emission of 
the last cycles of inspiral, followed by the merger and 
ringdown phases. Furthermore, by properly combining 
PN predictions and NR results, it is possible to describe 
analytically and/or numerically with high accuracy, the 
full gravitational- wave signal [5-8] . 

Black holes in binary systems can carry spin, and when 
spins are not aligned with the orbital angular momen- 
tum, spins induce precession of the orbital plane (see e.g. 
Rcf. [9]). This adds substantial complexity to the grav- 
itational waveforms, making it indispensable to include 
spin effects in templates used for the search. Moreover, 
as found long time ago [10-18], gravitational waves con- 
tain tail effects which are due to the back-scattering of 
linear waves in the curved space-time geometry around 
the source (and primarily generated by its mass). This 
causes the gravitational- wave signal to depend on the en- 
tire history of the binary system. 

In this paper we improve the knowledge and accuracy 
of the two-body inspiraling dynamics and gravitational- 
wave signal by computing the spin-orbit (SO) terms in- 
duced by tail effects. This is the continuation of our 
previous work on spins [19, 20] where we obtained the 
next-to-leading 2.5PN SO contributions in the equations 
of motion and gravitational-wave energy flux. Here, we 
derive those SO terms at 3PN order in the gravitational- 
wave energy flux, where they are entirely due to tails. 
Furthermore we obtain the SO terms at 2.5PN and 3PN 
orders in the wave polarizations that are specifically due 
to tails, leaving aside other SO terms at these orders that 
come from instantaneous (non-tail) linear contributions 
and which will not be computed here. 

We obtain the energy flux in two independent ways, 
first directly using the radiative multipole moments, and 
second by differentiating and squaring the gravitational- 
wave polarizations. To compute the SO tail effects in 
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the wave polarizations we solve the two-body dynam- 
ics taking into account spin precession. Assuming quasi- 
circular adiabatic inspiral, we also compute the 3PN SO 
terms induced by tails in the gravitational-wave phas- 
ing. These results can be used to improve the accuracy 
of inspiraling templates, to carry out comparison with 
numerical-relativity predictions, and to improve the ac- 
curacy of cffcctivc-onc-body and phenomenological tem- 
plates [5-8]. 

As an important check of our results we obtain the 3PN 
SO tail terms in the energy flux in the test-particle limit, 
and find perfect agreement with earlier PN computations 
based on black-hole perturbation induced by the motion 
of a test particle around a massive black hole [21]. 

This paper is organized as follows. In Sec. II we re- 
view the post-Newtonian multipole moment formalism 
and discuss relevant properties of tails. In Sec. Ill, we de- 
scribe how spin effects are included in the PN formalism 
and derive the binary's evolution equations when black 
holes carry spins. In Sec. IV we obtain the time evolution 
of the moving triad and solve the precessing dynamics at 
the relevant PN order. In Sec. V we compute the 2.5PN 
and 3PN SO tail effects in the gravitational waveform 
and polarizations. Restricting ourselves to quasi-circular 
adiabatic inspiral, we derive in Sec. VI the 3PN SO tail 
effects in the energy flux and in the gravitational phas- 
ing. Section VII contains our main conclusions. In Ap- 
pendix A we check that the equations of motion do not 
contain physical SO effects at 3PN order. In Appendix B 
we investigate the secular evolution of the spins, and fi- 
nally in Appendix C we give the explicit gravitational- 
wave polarizations. 



B. Notation 

In this paper we adopt the following notations. G is 
the Newton constant and c is the speed of light. As 
usual we refer to nPN as the post-Newtonian terms with 
formal order 0{c~^") relative to the Newtonian accel- 
eration in the equations of motion, or to the lowest- 
order quadrupole-moment formalism for the radiation 
field. Greek indices are space-time indices, and Latin are 
space indices. The quantity L = ii ■ ■ ■ ii denotes a multi- 
index composed of £ multipolar spatial indices ii, - ■ ■ ,ii 
ranging from 1 to 3. Similarly i — 1 = ii • • • ie-i and 
kL — 2 = kii ■ ■ ■ Nl ~ Ni^ ■ ■ ■ Ni^ is the product of 
£ spatial vectors Ni (similarly for xl = Xi^ ■ ■ ■ Xi^). We 
indicate with Ol = di^ ■ ■ ■ di^ and dkL-2 = dkdi^ ■ ■ ■ di^_^ 
the product of partial derivatives di = d/dx^. In the 
case of summed- up (dummy) multi- indices L, we do 
not write the £ summations from 1 to 3 over their in- 
dices. Furthermore, the symmetric-trace-free (STF) pro- 
jection is indicated using brackets, T^^^ = STF[Ti]; thus 
Ul = U(^L) and Vl = V(^l) for STF moments. For in- 
stance we write xi^iVj^ = \{xiVj + XjVi) — \5ijX ■ v. The 
transverse-traceless (TT) projection operator is denoted 



projector orthogonal to the unit direction AT = X /R of 
a radiative coordinate system = (cT, X). The quan- 
tity Eijk is the Levi-Civita antisymmetric symbol such 
that £i23 = 1. Finally, we denote time derivatives with 
a superscript (n), and we indicate the symmetrization 
operation with round parentheses. 



II. WAVE GENERATION FORMALISM 

The gravitational waveform hj^^ , generated by an iso- 
lated source described by a stress-energy tensor T''" with 
compact support, and propagating in the asymptotic re- 
gions of the source, is the TT projection of the metric 
deviation at the leading-order 1 /R in the distance to the 
source. It is parametrized by STF mass-type moments 
Ul and current- type ones Vl, which constitute the ob- 
servables of the waveform at infinity from the source and 
are called the radiative moments [22]. The general ex- 
pression of the TT waveform, in a suitable radiative co- 
ordinate system X'^ = {cT,X), reads, when neglecting 
terms of the order of 1/R^ or higher. 
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Here the radiative moments Ul and Vl are functions of 
the retarded time = T — R/c in the radiative coordi- 
nate system (we denote R= \X\). The integer £ refers to 
the multipolar order, and N = X / R = [Ni) is the unit 
vector pointing from the source to the far away detector. 
The TT projection operator Vj^^i and other notations 
are defined in Sec. IB. With P = {Pi) and Q = {Qi) 
denoting two unit polarization vectors, orthogonal and 
transverse to the direction of propagation TV. the two 
"plus" and "cross" polarization states of the waveform 
read as 



PiPj QiQj 7 TT 



, _ PjQj + PjQi , TT 



(2.2a) 
(2.2b) 



Our convention for the choice of the polarization vectors 
P and Q in the case of binary systems will be specified in 
Fig. 2. Plugging Eq. (2.1) into the standard expression 
for the gravitational-wave energy flux we get [22] 
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A. Expression of the radiative moments 

In the multipolar-post-Minkowskian formalism [IS- 
IS], the radiative moments are expressed in terms of two 
other sets of moments, referred to as the "canonical" mo- 
ments Ml , Sl , and which are relevant to the description 
of the source's near zone. The relation between the ra- 
diative moments J7l, Vl and the canonical ones Mi, Sl 



encodes all the non-linearities in the wave propagation 
between the source and the detector [15]. Those rela- 
tions may be re-expanded in a PN way and arc then seen 
to contain, at the leading 1.5PN order, the contribution 
of the so-called gravitational- wave tails, due to backscat- 
tering of linear waves onto the space-time curvature as- 
sociated with the total mass of the source itself. The 
explicit expressions at 1.5PN order arc [15, 23] 
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(2.4a) 
(2.4b) 



where AI is the Arnowitt-Descr-Misncr (ADM) mass as- 
sociated with the source. It also coincides with the 
mass monopole moment. The remainders C'(c~^),ion-taii 
in Eqs. (2.4) denote some correction terms which are at 
least of order 1.5PN and are instantaneous or contain the 
non-linear memory effect which will not be considered in 
the present computation. Here ni and tti denote some 
numerical constants given by [23] 
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(2.5a) 
(2.5b) 



The constant tq in Eqs. (2.4) is a freely specifiable time 
scale entering the relation between the radiative time Tr 
and the corresponding retarded time in harmonic coor- 
dinates. 

The canonical moments Ml and Sl arc themselves 
linked to six sets of multipole moments characterizing 
the source, collectively called the "source" moments and 
denoted /l, Jl, Wl, Xl, Yl, Zl- The point is that those 
source moments are known as explicit integrals extending 
over the pseudo-energy-momentum tensor of the matter 
fields and the gravitational field of the source [23, 24]. In 
the following we shall essentially need II and Jl which 
represent the main mass and current moments of the 
source. The other moments Wl, Xl, Yl and Zl play 
a little role because they are associated with a possible 
gauge transformation performed at linear order. It turns 
out that the difference between the canonical moments 
Ml-, Sl and the source moments II-, Jl arises only at 
the 2.5PN order: 
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(2.6a) 
(2.6b) 



For the present application it will be sufficient to con- 
sider the source mass moments II at IPN order and the 
current ones Jl at Newtonian order (see the discussion 
in Sec. IIIB). These are given by [25] 
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(2.7a) 
(2.7b) 



The other moments we shall need are the mass monopole 
moment M and the monopole of the moment Wl , which 
are given by 



M = / 



O 



1 



d^xx'a, + 



(2.8a) 
(2.8b) 



The mass, current and tensor densities a, <7i, cry 
Eqs. (2.7)~(2.8), are defined as (where T" = 5ijT'i) 



r 



Oi 



(2.9) 



We recall that, e.g., cr'-"-' in Eqs. (2.7) means taking n- 
time derivatives. 

The spin parts of the source moments in Eqs. (2.7) 
will come from the model we adopt for the stress-energy 
tensor T^"^ appropriate to spinning compact binaries (see 
details in Sec. III). Importantly, we notice that for the 
accuracy required by our calculation of the spin effects 
due to tails all integrands in Eqs. (2.7) have compact 
support. This is in contrast with the spin effects at 2.5PN 
order which necessitate non-compact supported higher- 
order terms in the source moments [20]. (We also find 
that the second term in II can be ignored in the present 
application to spins.) 
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B. Computing the tail integrals 

The tail integrals in Eqs. (2.4) extend over the en- 
tire past of the evolving source and it is a priori a non 
trivial task to compute them. Here we recall, based on 
Refs. [15, 16], that the tails are actually very weakly sen- 
sitive (in a post-Newtonian sense) to the past history of 
the source, and can essentially be computed using the 
current dynamics, i.e. at current time Tr, of the source. 

We have to compute, e.g., the integral appearing in the 
radiative mass multipole moment (2.4a), in which we can 
replace, following (2.6a), the canonical moment Ml by 
the source moment I^. Thus, 



J — OO 



r\t) In 



2r^ 



(2.10) 



where we pose Tq — tqc^'^'^. 

Let us introduce a constant time interval T to split 
the integral (2.10) into some contribution coming from 
the "recent past" , and extending from the current time 
Tr to Tr — T, and the remaining contribution called the 
"remote past", from Tr — T to — oo in the past. The 
recent past can be thought of as corresponding to the 
most recent orbital period of a compact binary system, 
while the remote past will include the details (eventually 
unknown) of the formation and early past evolution of 
the compact binary. However we shall prove that our 
result is independent of the chosen time scale T. 

To control the convergence of the tail integral (2.10) 
in the past we make a physical assumption regarding the 
behavior of the multipole moment Ihit) when t — > —oo. 
We assume that at very early times the source was formed 
from a bunch of freely falling particles initially moving on 
some hyperbolic-like orbits, and forming at a later time a 
gravitationally bound system by emission of gravitational 
radiation. The gravitational motion of initially free parti- 
cles isgivenby a;'(t) = VH+Wnn{-t)+X'' + o{l), where 
and denote constant vectors, and = GmV^/V^ 
(see Ref. [26] for a proof; to simplify we consider the rel- 
ative motion of two particles with total mass m). Here 
the Landau remainder o-symbol satisfies 9"o(l)/9t" — 
o(l/i"). From that physical assumption we find that the 
multipole moment behaves when t — > —oo like 



hit) = ALf + BLt'^'\ni-t) + CLt 



-o{t' 



(2.11) 



where A^, Bl and Cl are constant tensors. The time 
derivatives of the moment appearing in Eq. (2.10) are 
therefore dominantly like 



{t)^DLt-'' + o{t~''), 



(2.12) 



which ensures that the integral (2.10) is convergent. 

Next we integrate the "remote-past" integral (from 
Tr — T to —oo) by parts and make use of our assumption 
(2.11)-(2.12) to arrive at (posing t = Tr - Tx) 
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At this stage it is convenient to perform a Fourier de- 
composition of the multipole moment, i.e. 



, Zn 



(2.14) 



(The Fourier coefficients satisfy /£(f^) = Il{—^) since 
the moment is real.) Inserting (2.14) into (2.13) we ob- 
tain a closed-form result in the Fourier domain thanks to 
the mathematical formula [27] 



A/ dxlnxe'^^ + i 
/o Jl 



-f oo 



dx 
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^|s(A)-i(lnlAl+7E) , (2.15) 



where A = flT, with s(A) and ]A] denoting the sign and 
the absolute value, and 7e being the Euler constant. Fi- 
nally the result reads 



ULiTR) = if °^^i-in)'+'i^(n)e-'^^- 

J-oo 27r 

X |s(l])+i(ln(2ir!lr^)+7E)] . (2.16) 



We observe that the arbitrary time scale T has cancelled 
from this result. 

Later we shall apply this result to the computation 
of the waveform and energy flux of a spinning compact 
binary system. A priori^ since the tail integral (2.10) 
depends on all the past history of the binary (with the 
current binary's dynamics being the result of its long 
evolution by gravitational radiation emission) , we expect 
that the binary's continuous spectrum of frequencies 51 
should contain all orbital frequencies at any epoch in the 
past, say uj{t) with t < Tr, besides the current orbital 
frequency lo{Tr). However, it has been shown in the Ap- 
pendix of Ref. [16] that one can actually compute the 
tail integral by considering only the current frequency 
ui{Tr). Indeed the error made by this procedure is small 
in a post-Newtonian sense, being of the order of 0(^ In^), 
where ^ = lu/lo^ denotes the adiabatic parameter associ- 
ated with the gravitational radiation emission, and eval- 
uated at the current time Tr. In a PN expansion we have 
^{Tr) — 0(l/c^) so the error made by replacing the past 
dynamics by the current one is of the order of O (In c/c^) 
and can be neglected. The proof given in Ref. [16] is 
based on a simple model of binary evolution in the past, 
where an always circular orbit is decaying by radiation 
following the lowest order quadrupole formula, and spins 
are neglected. In this paper we shall assume that this 
result remains valid for spinning binaries. 
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III. APPLICATIONS TO SPINNING BINARIES 
A. Spin vectors for point-like objects 

Following our previous work [19, 20] we base our cal- 
culations on the model of point-partieles with spins [28- 
40] . The stress-energy tensor T''" of a system of spinning 
particles is the sum of a monopolar piece, made of Dirac 
delta-functions, plus the dipolar or spin piece, made of 
gradients of delta-functions: 
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, , [p u)5^'^\x-yA) 
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\/-9A 

p(m ^.i') S'-' ^^x - va) 
V--9A 



(3.1) 



where (S'^^ is the four-dimensional Dirac function, is 
the field point, y'^ is the world-line of particle A, Uj^ = 
dy^/icdTA) is the four- velocity, such that g^^u'^u'^^ = — 1 
where g"^^j = g^uiVA) denotes the metric at the parti- 
cle's location, is the linear momentum of the particle, 
and S'^ denotes its antisymmetric spin angular momen- 
tum. Our notation and conventions are the same as in 
Refs. [19, 20] which provide more details, except that here 
we shall denote using an overbar (i.e. Sj^) the original 
spin variable used in [19, 20]. Note that with our con- 
vention the spin variable has the dimension of an angular 
momentum times c. 

In order to fix unphysical degrees of freedom associated 
with an arbitrariness in the definition of S^'^ in the case 
of point particles (and associated with the freedom in the 
choice for the location of the center-of-mass worldline of 
extended bodies) , we adopt the covariant supplementary 
spin condition also called Tulezyjew condition [33, 34]: 



CM!' A 
'^A Pf 



0. 



(3.2) 



which allows the natural definition of the spin four- vector 
in such a way that 



\J-gA 



rriAC 



(3.3) 



where s^^^" is the four-dimensional antisymmetric Levi- 
Civita symbol (such that gOi^s — i)_ Por the spin vector 
itself, we choose a four-vector that is purely spatial 
in the particle's instantaneous rest frame, which means 
that in any frame 



(3.4) 



This choice is also adopted in Refs. [41-44]. As a eon- 
sequence of the condition (3.2), we can check that the 
mass defined by m^c^ = —p^P^ is constant along the 
trajectories, i.e. dm^/dr^ = 0. 

Important simplifications occur in the case of SO in- 
teractions, which are linear in the spins. Neglecting 
quadratic (spin-spin) interactions, the linear momentum 



is simply linked to the four velocity as = m^cu^, 
so the supplementary spin condition (3.2) reduces to 
S'a = 0, and the equation of evolution of the spins is 
given by 



-^ = 0, 
dTA 



(3.5) 



which means that the spin is parallely transported along 
the particle's trajectory. 

Following [19, 20] we adopt in a first stage as the vector 
spin variable the contravariant components of the vector 
S\, which are obtained by raising the index on Sf by 
means of the spatial metric 7Y, denoting the inverse of 
the covariant spatial metric = evaluated at point 

A (i.e. such that "f^ltj ~ ^j)- Hence our initial spin 
variable is 



lA 



(3.6) 



This definition of the spin vector Sa = {S\) agrees with 
the choice already made in Refs. [43, 44]. 

At the leading SO approximation, the contravariant 
spin variables Sa defined by Eq. (3.6) coincide with the 
spin variables with constant magnitude broadly used in 
the literature (see, e.g., Ref. [42]). At the next-to- leading 
order, the variables Sa differ from constant magnitude 
spins and their relationship has been worked out at 2PN 
order in Eq. (7.4) of Ref. [20]. In the present paper we 
shall denote the constant magnitude spins by Sa (al- 
though they were denoted S^ in Refs. [19, 20]). We know 
that it is actually better when presenting final results to 
switch to the constant magnitude spins Sa since they 
have a simpler precession equation (and turn out to be 
secularly conserved, i.e., over a radiation-reaction time 
scale; see Ref. [45] and Appendix B below). 

For two bodies {A — 1,2) the relationship between the 
constant magnitude spins and the original spin variables 
up to IPN order is: 



Si^Si- 



1 



--(viSi)vi^ Si 

2 ri2 



-0 



(3.7) 



together with the relation for the other particle obtained 
by exchanging all particle labels. We denote by va = 
dyA/dt the coordinate velocity of the particle A (with 
mass itia) and by ri2 = \yi — 2/2! the relative distance. 
See Ref. [20] for more accurate formulas extending (3.7) 
to 2PN order. In the present paper we shall consistently 
work only with the constant magnitude spins Sa- 
In the case of binary systems it is convenient to pose 



where Xi 



S^Si^ 
S2 

mi I'm and A2 



S2_ 

" Xi ' 

: TO2/TO (with m 



(3.8a) 
(3.8b) 



mi 



TO2). In addition we find it useful to occasionally use the 
dimensionless spin variables 

S S 



s = 



G m^ 



(3.9) 
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B. Multipole moments with spin-orbit effects 

The matter-source densities (2.9) depend on the com- 
ponents of the stress-energy tensor. At the leading PN 
order, the spin contribution therein (indicated by the 
subscript S) reduce to 



z JIJ 



fcfjS'^9fc(5i + lo2 + C'(^) , (3.10a) 

(3.10b) 

-^ekl(^v{^ S'ldi6i + l^2 + o(^^) , (3.10c) 



^3 



s 2c \c-^ 



where Si{x, t) = 6[x — yi{t)] means the three-dimensional 
Dirac delta- function evaluated on the particle 1, and 
1 <-> 2 means the same quantity but corresponding to 
the particle 2. 

In Ref. [20] the SO terms have been computed in the 
source mass quadrupole moment up to next-to-leading 
2.5PN order and the source current quadrupole moment 
Jij up to next-to-lcading 1.5PN order. All the other 
source moments were computed at the leading SO order. 
Those results are sufficient for our purpose. Actually, to 
compute the specific contributions of tails we need only 
the moments at leading SO order, given for general £ by 



II 

s 



Jl 

s 



21 



c3(£+l) 



(3.11a) 



L-2) 



1 O 2 



oil 



2c 



1^2 



c-^ 
(3.11b) 



Because the leading SO terms scale as C'(l/c'^) in the 
mass source moments, and as 0(l/c) in the current 
source moments, the number of non-linear terms needed 
in the radiative moments [Eqs. (5.1) below] is small. We 
refer to Sec. V of [20] for higher-order expressions of SO 
contributions of the source quadrupole moments. 



C. Equations of motion with spin-orbit effects 

Here we investigate the case where the binary's orbit is 
nearly circular, i.e., whose radius is constant apart from 
small perturbations induced by the spins (as usual we 
neglect the gravitational radiation damping at 2.5PN or- 
der). We denote by a; = yi — y2 the relative position of 
the particles (and v ~ dx/dt). Following Ref. [42] we in- 
troduce an orthonormal moving triad {n. A, £} defined by 
n = x/r a.s before, £ = Xn/I^^nI where Ln = rnvx x v 
with V = X1X2 denotes the Newtonian orbital angu- 
lar momentum and v the symmetric mass ratio, and 
A = £xn. Those vectors are represented on Fig. 1, which 
shows the geometry of the system. The orbital frequency 
w is defined for general, not necessarily circular orbits, by 
V = fn + ruX where r represents the derivative of r with 




FIG. 1: We show (i) the source frame defined by the orthonor- 
mal basis {x,y,z), (ii) the instantaneous orbital plane which 
is described by the orthonormal basis {xe, ye,£), (iii) tire mov- 
ing triad (n, X,£), and (iv) the direction of the total angular 
momentum J (agreeing by definition with the z-direction) . 
Dashed lines show projections into the x~y plane. 



respect to the coordinate time t. It is also equal to the 
scalar product of n and v which we denote as (nv) = r. 
The components of the acceleration a = dv /dt along the 
basis {n, A, are then given by 



n • a 
A • a 

£ ■ a 



ruj , 
+ 2rw , 



d£ 



(3.12a) 
(3.12b) 

(3.12c) 



We project out the spins on this orthonormal basis, defin- 
ing S = SnTi + S\\ + Se£ and similarly for S. Next 
we impose the restriction to quasi-circular processing or- 
bits which is defined by the conditions f = = r so 



that v"^ 



2 2 
r uj 



(neglecting radiation reaction damping 



terms). In this way we find [19] that the equations of the 
relative motion in the frame of the center-of-mass are 



dv 



UJ n - 



(3.13) 



There is no component of the acceleration along A. Com- 
paring with Eqs. (3.12) in the case of circular orbits, we 
see that uj is indeed the orbital frequency, while what 
we call the "precessional frequency" Wproc = A • d£/dt is 
proportional to the variation of £ in the direction of the 
velocity v = rujX. We know that is given by 



G? 



1 + 7 (-3 + i^) + 7^/2 {-5se - Sdai) 



(3.14) 



where we denote S = Xi — X2 and S( = {s£) = s ■ £, 
where the spin variables are defined by Eq. (3.9). The PN 
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parameter is 7 = Gm/{rc^) and we have included only 
the IPN non-spin term and the leading SO correction at 
1.5PN order. On the other hand, we get [19] 



3/2 



7s„ + iScTj 



O 



(3.15) 



where s„ = {sn) = s ■ n. At the leading 1.5PN SO 
order the orbital frequency (3.14), as well as Wprcc, remain 
unchanged if we were to substitute some other variables 
to the spins S, S. However, when working at a higher 
FN approximation, it is more convenient to use the spin 
variables S, S, built from the constant magnitude spins. 
The main advantage of the spins Sa is that they satisfy 
the usual-looking precession equations 



dSA 
dt 



flA X S/ 



(3.16) 



showing that the spins precess around the direction of 
Qa, and at the rate ft a = \^a\- The equation (3.16) 
could in principle be extended to any FN order (at the 
linear SO level). The precession's angular- frequency vec- 
tors flA have been computed up to the 2FN order for 
circular orbits in Ref. [20]. Here, we shall only need the 
IFN leading order: 



£ + 



1 



(3.17) 



To obtain we simply have to change S into —S. Both 
precession frequencies are constant in magnitude and in- 
dependent of the spins in the 1.5FN dynamics. 

The equations of motion (3.13) and the precession 
equations (3.16) together leave invariant the total an- 
gular momentum. 



c 



(3.18) 



where L denotes the orbital angular momentum. For 
future reference we give the components of L along the 
triad basis at IFN order for non-spin effects and at the 
leading 1.5FN order for spin ones [19, 20]: 



c 



.-1/2 



1 + 






(2^ 











Sri 



1 



2 2 

-3Sx - 



(3.19a) 
(3.19b) 
(3.19c) 



Note that the components L„ and Lx are due to spin 
effects arising at order 0{c-^). See Eq. (7.10) of Ref. [20]. 



IV. EVOLUTION OF THE TRIAD {n,\,e.} 

Using Eq. (3.13) the time derivatives of the three mov- 
ing triad vectors {n. A, £.} can be expressed with respect 
to that triad basis as 



dn 
dF 



(4.1a) 




FIG. 2: Similar as Fig. 1 but with the direction of the source 
N indicated together with a choice of convention for the two 
polarization vectors P and Q. 



dX 

M 

dt 



■^prcc ' 



(4.1b) 
(4.1c) 



Equivalently, introducing the orbital precession vector 
— UJ £ and spin precession vector Wpj-oc = — '^proc'T-, 
these equations can be written as 



dn 

lit 
dX 

'dt 
d£ 

di 



a; X n , 

prcc) X A , 

x£. 



prcc 



(4.2a) 
(4.2b) 
(4.2c) 



We recall that the spin precession frequency is given by 
Eq. (3.15) or equivalently 



7s„ + S^df 



O 



(4.3) 



where we pose x = {Gmto / c^)"^/^ which defines a gauge- 
invariant FN parameter agreeing with 7 at leading order. 

We shall now solve the evolution equations for the mov- 
ing triad {n, A, £} at the 1.5FN order in the conservative 
dynamics (i.e., neglecting radiation reaction). We recall 
that the spin variables we use in this calculation are those 
with constant magnitude, denoted here Sa- 

It is convenient to introduce a fixed (inertial) orthonor- 
mal basis {x,y,z}; see Figs. 1 and 2. For a given value 
of the total angular momentum J (a constant vector), 
and of the direction N = X/R oi the detector as seen 
from the source, a canonical choice of the basis vectors is 
as follows: (i) z is defined to be the normalized value of 
namely J/ J; (ii) y is orthogonal to the plane spanned 
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by TV and z and points to the direction that corresponds 
to the positive orientation of the acute angle (2;, TV), i.e. 
y=:2;xTV/|2;xTV|; (iii) x completes the triad. We see 
that a;, z and TV are coplanar by construction. Then, 
we introduce the standard spherical coordinates with the 
inclination angle measured from the zenith direction z 
and the azimuthal angle measured from x. The spher- 
ical coordinates of TV and I are denoted as (6*, <~p) and 
(a, i) respectively, and since TV lies in the same plane as 
X and 2, we have = (see Fig. 2). Since t is the angle 
between the total and orbital angular momenta, we have 



sm L 



I J X £| 
J 



(4.4) 



The angles (a, l) are referred to as the precession angles. 

We now derive the time evolution of our triad vectors 
from that of the precession angles (a, t), and of an ap- 
propriate phase $ that specifies the position of n with 
respect to some reference direction. Following Ref. [46], 
we introduce the unit vectors 



J X £ 
|Jx£| ' 



X Xi , 



(4.5) 



such that {xi,yi,£.} is an orthonormal basis. The phase 
angle $ is defined by (see Fig. 1): 



{xe,n) 



(4.6) 



The rotation takes place in the instantaneous orbital 
plane spanned by n and A, and we have 



n ~ cos ^ xi + sin $ yi , 
A = — sin <t> ajf + cos $ ye , 



from which we deduce 

Xi ■ in 



iA = 



Combining (4.8) with (4.4) we also get 



sm L e 



J 



(4.7a) 
(4.7b) 



(4.8) 



(4.9) 



By identifying the right-hand sides of Eqs. (4.1) or 
(4.2) with the time-derivatives of the identities (4.7) we 
obtain the following system of equations for the varia- 
tions of a, i and $, equivalent to the system (4.1), 



da 

'dt 

di 

di 
d$ 
'dt 



sin$ 



^prcc 



^prcc 



sm L 
COS$ , 

sin$ 



^prcc 



tan ( 



(4.10a) 
(4.10b) 
(4.10c) 



On the other hand, using the total angular momentum 
(3.18) together with the components of the orbital angu- 
lar momentum given by Eqs. (3.19) — notably the fact 



that Ln and L\ are due to SO terms dominantly of order 
0{c~^), we deduce that sint is a small quantity of order 
0{l/c). From this fact, we conclude by direct integration 
of the sum of Eq. (4.10a) and Eq. (4.10c) that 



$ + a = (/) + C'(^) 



in which we have defined the "carrier" phase as 



u dt = uj[t — to) + : 



(4.11) 



(4.12) 



with 00 the value of the carrier phase at some arbitrary 
initial time to. We recall that the orbital frequency (3.14) 
is constant in first approximation for circular motion. 

The combination $ -I- a being known by Eq. (4.11), 
we can further express the precession angles i and a in 
first approximation in terms of the components Sn and 
S\ of the total spin S ^ Si + 82- From (4.4) we find 
(discarding non-linear spin contributions) 



SI 



O 



(4.13) 



where we recall that Ln = mvr^uj denotes the Newtonian 
orbital angular momentum. On the other hand, using 
also Eq. (4.9) and the relation (4.11) we obtain at leading 
order 



S\~\Sn 



SI 



o 



.14) 



[See also the more precise equations (4.21)-(4.22).] 

It remains now to obtain the explicit time variation of 
the components of the individual spins 5*^, and Sf. 
Using (4.13) and (4.14) [and also (4.11)] we shall then be 
able to obtain the explicit time variation of the precession 
angles and phase. Combining (3.16) and (4.1) we obtain 
the precession equations for the three unknowns , 
and sf in the form of the following first-order system 
(valid at any FN approximation) 



dsl 

dt 

dSl 

dt 
dSf 

dt 



nA)S^ 



^A)Sf 



(4.15a) 
(4.15b) 
(4.15e) 



where $7^ is the norm of the precession vector of the spin 
A as given by (3.17), and the precession frequency Wpioc 
is explicitly given by (4.3). Actually the terms involving 
Wpioc in the right-hand sides of (4.15) can be neglected 
because they are quadratic in the spins. Thus, staying 
at the linear SO level, we find that the equations (4.15) 
can be decoupled and integrated as 



Sj_ cos IpA 



St^-Sfsmi^A 



(4.16a) 
(4.16b) 
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Si 



(4.16c) 



Here and S-f^ denote two constants for each spins A, 
and agree with the projections (which are constant at the 
Hnear SO level) of the spins perpendicular and parallel 
to the direction of the orbital angular momentum £.. The 
phase of each of the spins is given by 



V^A = {UJ - nA){t - to) + 1p°A 



A > 



(4.17) 



where i/jJ^ is the constant initial phase at time t^. 

With those results we obtain an explicit solution for 
the precession angles by substituting Eqs. (4.16) into the 
results (4.13) and (4.14). We find that i(t) is given at 
the 0.5PN level by 



4/2 



sm L = 



(4.18) 



O 



where we recall that = S^/{Gm'^). Knowing i(i) we 
deduce a{t) from 



r.1/2 



sm i e 



s2 e-"^= 



O 



(4.19) 

The difference of spin phases 4'i2 = '0i~'02 readily follows 
from Eq. (4.17) and Eq. (3.17) at IPN order as 

^-12 = <2 + ^^a-<5(t-to)+o(^) . (4.20) 

Moreover, Eq. (4.9) can be written more explicitly at the 
1.5PN level as 



sm L e 



Li Vc4 



(4.21) 



where J_|. = J„ + i Ja is given at the 1.5PN order by 

S] 



J4 



-1-01 



,13 6 



(4.22) 



S ly 35 



and the IPN orbital angular momentum is known 
from Eq. (3.19a). 

As a check of the previous solution we observe that if 
we take the time derivative of Eq. (4.4), then evaluate 
the total angular momentum J given by (3.18) together 
with the components of the orbital angular momentum 
L provided in (3.19), and use the solution (4.15) for the 
evolution of the spin components, we obtain 



di 



5, 



(4.23) 



which is consistent with (4.10b) once (4.8) is employed. 

Finally we express the triad vectors n(t), \{t) and £{t) 
in terms of the precession angles and the carrier phase, 
and in terms of the initial triad and angles at the initial 
instant tg, modulo terms of order 0{c~'^). To do this 
we notice that the triad {n. A, £} at time t is obtained 
from the inertial triad {x, y, z} by the rotation associ- 
ated with the three Euler angles a, l and $. Similarly 
the initial triad {noj-^Oj-^o} at time is obtained by 
the rotation associated with ao, to and $o- So, com- 
bining those two rotations we readily obtain {n. A, in 
terms of {ng, Ao,£o}- Using Eq. (4.11) to eliminate the 
phase $ in favor of the carrier phase (j) — this introduces 
small remainder terms 0{c~'^) — and neglecting all terms 
quadratic in the spins, we get 



n = cos 



4>o)no + sin((/) — 0o)Ao + ^sin t sin(0 — a) — sinio sin(0 — ao)^-^o + c(^— ^ , (4.24a) 
A = — sin(0 — (/)o)no + cos((/) — (j)o)\o + |^sin tcos(0 — a) — sinio cos(0 — ao)^-^o + o(^—^ , (4.24b) 
£ = £0 + (^— sintsin((/)o — a) -I- sinio sin((/)o — ao)^no + (^— sinicos((/)o — a) + siuio cos((/)o — ao)^ Ao + O^— ^ . 



(4.24c) 



Since we have found in Eqs. (4.18)-(4.19) an explicit so- 
lution for the precession angles i(t) and a(t), the time 
dependence of n, A and £. is completely known. Note 
that in practical computations it is often more conve- 
nient to work not with n and A but with the complex 
null vector m = {n + i A)/\/2 and its conjugate. 



V. COMPUTATION OF THE WAVEFORM 

Here we shall compute the SO terms coming from all 
non- linear (i.e. of formal order G^) contributions associ- 
ated with tails consistent with the 2.5PN and 3PN orders 
in the waveform. We shall need only to focus on the tails 
entering the mass and current quadrupolcs Uij and Vij 
(having £ = 2) and on the current octupole Vijk (^ = 3). 
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They indeed contain, when speciahzed to spinning com- 
pact binary systems, the SO contributions we are inter- 
ested in. The reason is that the leading SO terms start 
at the 0.5PN order 0{l/c) in the cmrent moments, but 
only at the 1.5PN order 0{l/c'^) in the mass moments 
[see Eqs. (3.11) above]. 

In addition to the tail integrals shown in Eq. (2.4), we 
shall also compute some terms of order at 2.5PN or 
1.5PN order, but which are non- hereditary, i.e., merely 
depend on the instantaneous retarded time Tr. Those 
corrections are given in full form by Eqs. (5.4)-(5.6) of 
Rcf. [47] , but here we shall need only, for the same reason 



as before, two terms involving the source current dipole 
moment or angular momentum J^; the other terms will 
not contribute to the SO effect at 3PN order. Further- 
more, we have to include some additional corrections of 
similar nature originating from the differences between 
the canonical and the source moments in Eqs. (2.6). 
Those arc given in the general case by Eqs. (5.9)-(5.10) 
of Ref . [47] , but we can check that only the contribution 
in the mass quadrupole Uij a priori matters here. 

The relevant contributions in the radiative moments, 
including only the terms needed for the applications be- 
low (see Ref. [47] for more complete expressions) read: 



sv,. 



ijk 



2Gm 
2Gm 
2Gm 



dt 



dt 



dt 



R 



2to 

Tn-t 

2ro 
Tn-t 

2to 
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+ 12. 

7' 
+ 6 
5 

+ 3 



r(4) 
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(2) 



Ji5)(.^ 2G (4) 



(5.1a) 
(5.1b) 
(5.1c) 



where we have replaced M with m = mi + 7712 which 
is valid at the dominant order. Moreover, we shall find 
that the terms in the right-hand side of Eq. (5.1a) which 
depend on the moment W vanish at the considered order. 
Indeed, by inserting the value of di from (3.10b) into the 
potential W defined by (2.8b) and integrating by part, 
we find that the result is zero. 

The corresponding gravitational-waveform, for which 
all three moments Uij, Vij and Vijk are important, is 
then given by 



TT 



-pTT 

c^R 

1 

~ 2? 



Nam £ab{k ^^)bm 



(5.2) 



The main task consists in evaluating the tail integrals 
(2.10) whose Fourier transforms have already been ob- 
tained in Eq. (2.16). The result (2.16) heavily relied on 
a physical assumption concerning the system in the re- 
mote past, namely that it was formed by freely falling 
incoming particles, see Eq. (2.11). 

We reviewed in Sec. II B that one can insert in the re- 
sult (2.16) the binary's current frequency spectrum, i.e. 
at time Tr, modulo small error terms of the order of the 
adiabatic parameter of the inspiral, or, more precisely, 
of negligible order 0{\nc/c'). In the case of spinning 
compact binaries this means that we have to include in 
the spectrum the current orbital frequency w = w(Tr), 
and also the precession frequencies fli = fli(Tfi) and 
= il2{Tu) of the two spins. This follows from the ex- 
plicit solution of the triad {n. A, £} and of the precession 



equations (see Sec. III). [Notice that the precession angles 
a and l always appear through the product sin t e'" given 
by equations such as (4.19).] Hence we can take for the 
Fourier components of the multipole moments 



^ n.ni ,719 



(5.3) 



where the frequency modes are some w„.„j_„2 = nuj 
ni rti + n2 ^2- The result (2.16) then becomes 



ULiTn)^ ^ iAr-"^(-ic.„..„„J^+ie-'- 

n,ni ,712 

(2| i,"2 l''"o) ^" 7E I 



(5.4) 



We recall from Eq. (3.17) that the precession frequencies 
ill and ^2 are small quantities of order IPN. This means 
in particular that because of the explicit factor uj^ni.n^ 
in Eq. (5.4) (which arises from taking the time derivatives 
of the multipole moment and integrating), the modes for 
which n = in tail integrals are very small, at least of 
order 4.5PN, and can be neglected. 

The SO terms in the radiative tails originate primarily 
from the spins present in the sources of the integrals, i.e. 
appropriate derivatives of multipole moments as shown 
in Eqs. (5.1). The SO contributions in the source mo- 
ments have been given in Eqs. (3.11). However there 
are other crucial SO terms which originate from the non- 
spin parts of the moments. They come from time differ- 
entiations of the triad vectors using the evolution equa- 
tions (4.1), which produce spin terms contained in the 
precession frequency Wpicc, and from time differentiations 
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of spins themselves via the precession equations (3.16). 
SO contributions may also be generated by the tail in- 
tegration itself due to the precession of the triad basis 
{n, A, according to the formula (4.24). On the other 
hand, other SO terms come from the reduction to circu- 
lar orbits when we eliminate the orbital separation r in 
favor of a function of the orbital frequency uj obtained 
from inverting the relation (3.14). Note that the two lat- 
ter corrections being of 1.5PN relative order, they cannot 
come from anywhere but the tail integral of the Newto- 
nian quadrupole moment. 

During the practical computation, we make explicit 



J 



the time dependence of the derivatives of multipolc mo- 
ments, computed in the center-of-mass frame as func- 
tions of the relative position, the relative velocity and 
both spins. For circular orbits, x and v depend only on 
r and cj, which are approximately constant on dynamical 
time-scales, and on the unit vectors n and A. Thus, the 
whole time dependence arises through that of n and A 
(and £ = n X A), and is provided by our explicit solu- 
tion (4.24), together with the precessing angles a{t) and 
L{t) given by Eqs. (4.18)-(4.19), or (4.21) and (4.22) with 
more precision. 



The complete results for the spin dependent parts of the radiative moments, in which we use the short- hand notation 
for spins (3.9) and where the basis vectors {n,X,£} are evaluated at the current time Tr, are then 
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35(5cr,0 iyTTn^'e^ - 2X^'P'^ (ln(4a;ro) 
i(29sA + n5cjx) [t: A(V^) + 2n(*F) (ln(4wTo) + 7e - 





+ 4(s£ 




- 2^n 






5V,,k = 





A<V^'> fln(2a;To) 



7e 



19^ 
2 



3(5(7, 



(5.5a) 
(5.5b) 
(5.5c) 



Insertion of the above quantities into Eq. (5.2) yields the non-linearly induced SO contributions at 2.5PN and 3PN 
orders in the waveform as 



5h. 



TT 



.5/2 



A'^(ln(2Toc.)+7E-^) -^n^ 



{N X cr)M,fc + N''eak^a' 



: (^/-(n^^ - A^^) + 2(/^ + 5a'') ((n^-^ - A^"^) (ln(4roc<j) + 7e - ^) H 
+ 4(s, + ^) (-4n*A^ (ln(4Toc^) + 7e - ^) + 7r(n*^ - A*^)) - ^s,n'\^ 



2(n'e (^sx + 35 ax) + X'f f + 3(5ct„ 
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+ i(71s„ + 355an)(nn'£^ - 2V£-'' (ln(4c^To) + 7e - ^ 
- i(29sA + 175ax) {nyf + 2n'f (ln(4cjro) + IE - ^ 



TT 



(5.6) 



for which we have conveniently introduced the rescaled waveform 5hJ^ defined by 



(5.7) 
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The two gravitational-wave polarizations /i+ and h X are 
given in Appendix C. We have checked that the test- 
particle limit 1/ — > of the —2 spin- weighted spherical 
modes resulting from the above waveform agrees with the 
results of Ref. [21] (given explicitly in Ref. [48]) based on 
black-hole perturbation theory. 



VI. ENERGY FLUX AND ORBITAL PHASING 

The case of the gravitational energy flux is simpler 
than for the waveform, notably because we need only 
the contributions from the mass and current quadrupole 
moments, i.e. 



G 



5 «j 



(1) 



45c2 *J 



.1) 



The 3PN SO effects in the energy flux have been com- 
puted in two different ways. In the first way, we com- 
pute the time derivative of the radiative moments Uij and 
Vij whose SO-tail contributions are given in Eqs. (5.5), 
and then square these radiative moments to get the flux 
(6.1). The second way is completely equivalent, but en- 
tirely done by hands. It consists of writing all the sep- 
arate pieces composing the energy flux (6.1), made of 
the coupling between some instantaneous moment (eval- 
uated at current instant Tr) times a hereditary tail in- 
tegral. The SO terms have to be included in either the 
instantaneous moment in front of the integral, or in the 
tail integral itself. This gives then several "direct" SO 
contributions coming from tails at relative 1.5PN order 
(for the mass quadrupole tail) or 0.5PN order (for the 



current quadrupole tail) which are then added together. 
In addition there is the crucial contribution due to the 
reduction to circular orbits of the standard (non-spin) 
tail integral at 1.5PN order, for which the relation be- 
tween the orbital separation r and the orbital frequency 
Lo [as given by the inverse of Eq. (3.14)] provides a sup- 
plementary SO term at relative 1.5PN order, which thus 
contributes in fine at the same 3PN level as the "direct" 
SO tail terms. 

Finally, we obtain the following net result for the SO 
tail contribution at 3PN order in the total energy flux: 



-167r Si 



3l7r 



5 ag 



(6.2) 



where we recall that S£ = s ■ t and og ~ cr ■ I, with 
the spin variables s and cr being defined by Eqs. (3.7)- 
(3.9). Let us remark that in the energy flux the 3PN 
SO term is entirely constituted by the SO tails we have 
obtained in (6.2). So the complete 3PN SO term in the 
flux is provided by Eq. (6.2). Contrary to the waveform 
computed in Sec. V, there are no other SO terms coming 
from linear source moments at that order. 

Because the energy flux and the resulting orbital phas- 
ing is so important for gravitational-wave observations, 
we shall now give the complete formula for the total flux, 
including all non-spin terms and all linear SO terms up 
to 3PN order (but neglecting non- linear SS interactions). 
However we shall not write the known non-spin 3.5PN 
terms in the flux (due to non-spin tails [17]) because some 
yet uncalculated SO effects should conjointly appear at 
that order. The 3PN energy flux, complete except for SS 
interactions, reads then 
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We are consistently using the constant-magnitude spins 
Sa that are related to the original variables Sa of 
Ref. [19, 20] by Eq. (3.7); see also Eqs. (7.4) of Ref. [20]. 
The non-spin terms are given, e.g., in Ref. [1]. We find 
perfect agreement in the pcrturbative limit — > with 



black hole perturbation calculations reported in Ref. [21]. 
On the other hand the total conservative energy E of the 
binary is not affected by the SO terms at the 3PN order 
(we check this point in Appendix A), hence we have 
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Following Ref. [20] we shall next use the standard en- 
ergy balance argument to deduce the evolution of the or- 
bital frequency even in the presence of spins. To this end 
we have to check that the constant-magnitude spins are 
secularly constant (i.e., constant over a long radiation- 
reaction time scale) up to the right level, 1.5PN order 
in the present case. In Ref. [20] we have referred to the 
work [45] for a proof that this is correct up to relative 
IPN order, i.e. considering radiation reaction effects up 
to 3.5PN order. In Appendix B below we extend the ar- 
gument of Ref. [45] to the relative 1.5PN order, which 
essentially means adding the tail-induced part of the ra- 
diation reaction at 4PN order. This check being done we 
can thus neglect (ds^/d<) and {dag/dt) in average over a 
radiation-reaction time scale. 

An alternative way to see this is to directly com- 
pute the variation of the projection of the spins along 
the Newtonian orbital angular momentum, i.e. Sf = 



Sa ■ using the precession equations (3.16) appropri- 
ate for constant-magnitude spins. We readily find that 
d5//dt = SA-[dl/dt+l-KVlA], which shows that AS^/dt 
is at least quadratic in the spins for circular orbits. This 
readily follows from the facts that i remains constant 
in the absence of spins, and that, as we have seen in 
Eq. (3.17), 0,A for circular orbits points in the direction 
of £ modulo spin corrections. Thus we have dS^ /dt = 
at the linear SO level (neglecting quadratic SS couphngs). 
The argument is in principle valid up to any PN order, 
but is restricted to circular orbits. 

The conclusion is that the constant-magnitude spin 
terms can be considered as constant when computing 
the averaged evolution {dE/dt) of the energy given by 
Eq. (6.4). Equating then (dE/dt) to —J", where T is 
given by (6.3), we obtain the secular variation of the fre- 
quency (lu) — denoted to for simplicity — as (neglecting 
SS contributions) 
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By integrating this out using standard PN rules for mul- 
tiplying, dividing and integrating PN expressions, we ob- 



tain the secular evolution of the carrier phase [defined by 
(/) = J ujdt; see Eq. (4.12)] as 



+ x 



1 

32^ 

-1/2 



r.-5/2 



„-3/2 



/ 15293365 



3715 55 
1 

12 

3085 



.1008 
27145 



V 1016064 1008 144 
/ 38645 554345 41745, 



In a; I --r— — tt 
V 1344 



2016 



448 



-IOtt 



65 



235 125 , 

Sf H 0(7/ 

6 8 



55 



16 



-Se 



15 



6(7p 



14 



- X 



1/2/12348611926451 
V 18776862720 

r 15737765635 
ly \ — — — — 



160 2 
2255 2 



1712 856 



21 



-7E 



21 



ln(16 x) - 



9407r 7457r 



-si 



12192768 ~^ 48 



76055 2 127825 

v 1/ 

6912 5184 J 



6 



■Sai 



(6.6) 



We recall that to the carrier phase we have also to add the 
precessional correction, arising from the changing orien- 
tation of the orbital plane. We have proved in Eq. (4.11) 
that at the IPN order the total phase $ is given by 
$ = — a + 0{c~'^). Thus the precessional correction 
is given by —a and is explicitly provided by the solution 
(4.18)-(4.19). Alternatively, the precessional correction 
can be computed numerically [9] . 



VII. CONCLUSION 

So far, the search for gravitational waves with LIGO 
and Virgo detectors has focused on non-spinning com- 
pact binaries [49-53], although in Ref. [54] single-spin 
templates were employed, for the first time, to search for 
inspiraling spinning compact objects. It is timely and 
necessary to develop more accurate templates which in- 
clude spin effects. Extrapolating results from the non- 
spinning case, we expect that, for maximally spinning 
objects, reasonably accurate templates would need to be 
computed at least through 3.5PN order. 

During the last years, motivated by the search for 
gravitational waves, SO effects have been computed in 
the two-body equations of motion through 3.5PN or- 
der [19, 41, 43, 44, 55-57] and energy flux through 2.5PN 
order [20, 42-45]. Moreover, SS effects have been cal- 
culated through 3PN order in the conservative dynam- 
ics [41, 42, 58-68] and multipole moments [69]. 

In this paper, building on our previous work [19, 20], 
we have improved the accuracy of the energy flux and 
gravitational waveform by computing SO terms induced 
by tail effects [10-18]. Those effects are due to the back- 
scattering of linear waves in the curved space-time geom- 
etry around the source. Using the multipolar PN formal- 
ism developed in Refs. [13-15, 23, 24], we have identified 
and computed the radiative multipole moments respon- 
sible of tail terms involving SO couplings. More specif- 
ically, we have computed those SO tail contributions to 
the energy flux at 3PN order and to the gravitational 
waveforms at 2.5PN and 3PN order. Those SO tails con- 
stitute the complete coefficient at 3PN order in the en- 
ergy flux. In particular we find that the energy flux is in 
complete agreement with the result of black-hole pertur- 
bations in the test-particle limit [21]. Our computation 
is restricted to quasi-circular inspiraling orbits, and uses 
the two-body precessional dynamics at 1.5PN order. 

The computation of SO tail terms in the waveform 
is summarized in Sec. V, and some building blocks and 
foundation for calculating tail effects in the PN formalism 
were reviewed in Sec. II. For the first time, we have 



computed tail terms when precession effects in the two- 
body dynamics are also included. The relevant results for 
the waveform are given in Eq. (5.6) and in Appendix C. 
The SO tail effects in the energy flux and phasing at 3PN 
order are given in Sec. VI, see in particular Eqs. (6.2)- 
(6.3) and (6.6). 

Considering the vigorous synergy which is currently 
taking place between analytical and numerical relativity 
for building faithful templates [5-8], we expect that the 
results developed in this paper will help the construc- 
tion of more accurate analytical templates describing the 
entire process of inspiral, merger and ringdown of black 
holes in presence of spins. 

In the near future we plan to complete the knowledge 
of SO effects in the gravitational waveform at 3PN order, 
by computing the non-tail (i.e., instantaneous) SO cou- 
plings at 2PN and 3PN orders, and the corresponding 
—2 spin-weighted spherical harmonics (or gravitational 
modes). This will constitute a step further with respect 
to Ref. [46] which computed SO effects in the gravita- 
tional modes through 1.5PN order. 
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Appendix A: 3PN spin terms in the equations of 
motion 

In this Appendix, we check that there are no SO terms 
at the 3PN order in the total conservative invariant en- 
ergy of the binary given by Eq. (6.4). Indeed, we find 
that the 3PN SO terms in the binary's equations of mo- 
tion (say, in harmonic coordinates) can be gauged away. 
The result is to be expected because we know that the 
first SO modification of the radiation reaction damping 
force arises at the 4PN order rather than 3PN [45]. 

We compute the near-zone PN metric by solving the 
Einstein field equations in harmonic coordinates for the 
stress-energy tensor (3.1). We find by direct PN iteration 
of the metric, parametrized by means of retarded poten- 
tials V, Vi, ■ ■ ■ (see Ref. [19] for more details), that the 
contribution of SO terms at 3PN order in this gauge is 
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where we indicate with cst(i) some irrelevant 0{c~^) con- 
stant term in space, where we keep the SO parts of the 
acceleration un-replaced, and where 1 2 refers to the 
same expression but for particle 2. 

The metric (Al) yields a 3PN contribution in the equa- 
tions of motion of spinning particles which can be calcu- 
lated from the Papapctrou [28, 29] equations of motion 
(see e.g. Sec. Ill in [19]). The result for the acceleration 
of particle 1 is 
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We observe that Sa\ is symmetric under the exchange of 
particles 1 and 2. A closer inspection reveals that it is in 
fact given by the second total time derivative of a certain 
vector, namely 



5a\ ~ 8a\ 



~dt2~ 



(A3) 



and we find that dX. = v^^(x x cr) in the notation 
of Sec. inc. This is precisely the effect of the gauge 
transformation associated with a shift of coordinates 
x'* = x'' + 6X\ We thus conclude that the 3PN SO terms 
in the equations of motion are pure gauge and cannot af- 
fect the binary's invariant energy (6.4). 



Appendix B: 4PN spin secular evolution 

Here we show that the constant-magnitude spins are 
secularly constant, i.e. constant over a long radiation- 
reaction time scale, up to the 4PN order corresponding 
to the 1.5PN relative order. In Ref. [45] this has already 
been proved up to IPN relative order; here we extend 
the argument to 1.5PN order. [In the main text after 
Eq. (6.4) we present an alternative argument valid at 
any PN order but restricted to circular orbits.] 

Following Ref. [45] we describe our source by a set of 
well-separated extended bodies A, supposed to be New- 
tonian in a first stage. We define the spin S\ of each 
of the bodies in the usual Newtonian way as an integral 
extending over the volume of body A, 



'ijk 



where denotes the Newtonian (baryonic) mass density, 
and x'^ is the Newtonian center of mass position of the 
body A. At Newtonian level the spin (Bl) agrees with the 
definition employed in the present paper. The "baryonic" 
spin (Bl) is only used for the purpose of this Appendix. 
The equation of evolution of the baryonic spin reads as 
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The spin precession equation follows from inserting 
into (B2) an explicit solution for the acceleration in terms 
of positions and velocities. The resulting equation is 
then simplified using some virial relations appropriate 
to the case where the compact body is "stationary" , see 
Ref. [45]. The secular evolution of the spin is then ob- 
tained by considering the radiation reaction piece of the 
acceleration in Eq. (B2). 

At the dominant 2.5PN level the radiation reaction 
acceleration inside an isolated body in harmonic coordi- 
nates is given by (see, e.g., Ref. [70]) 
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where lij is the source's STF quadrupole moment (at 



Newtonian order), and Uij 
tensor defined by 



is the Newtonian potential 



U^j{x,t)^G I d'x'p4x',t)^-^^^^^^^^^^. (B4) 



(We have Uu = U, the usual Newtonian scalar poten- 
tial.) It can be shown [45] that the only contribution at 
2.5PN order to the spin precession equation comes from 
the velocity-dependent part of Eq. (B3), i.e. 



2G 



^2.5PN - ^5 



(B5) 



The other pieces in the 2.5PN acceleration vanish when 
the size of the body tends to zero (compact-body hmit) 
and may be ignored. Using a virial relation [45] we read- 
ily obtain 
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Because the spin is constant in the lowest approximation 
the latter result is a total time-derivative: 
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which can be moved to the left-hand side and absorbed 
into a negligible redefinition of the spin variable at 2.5PN 
order. When specialized to two compact bodies the result 
(B7) becomes 
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[Note that the latter expression depends on the specific 
definition of the spin one is using, i.e. in the present case 
the Newtonian spin defined for extended bodies by (Bl); 
for the spin variable used in [19, 20] the expression is 
expected to be different, but still in the form of a total 
time derivative.] 

The results (B7) or (B8) show that there is no secular 
evolution for the spin at the dominant 2.5PN order (see 
also [71] and references therein for related discussions). 
Note that this conclusion actually applies to any spin 
variable at dominant order. However it has been shown 
in Ref. [45] that in the case of the constant-magnitude 
spin there is also no secular evolution of the spins at the 
next-to-leading 3.5PN order in the radiation reaction. At 
next-to-leading order this result [45] applies specifically 
to the constant-magnitude spins and uses the radiation 
reaction acceleration up to the 3.5PN order. 

We now extend the argument by including the higher- 
order 4PN correction term (i.e. 1.5PN radiation-reaction 
order) which is known to be due to gravitational wave 
tails [14, 72]. That extension is actually straightforward 
since it basically needs only the Newtonian result (B7). 
The reason is that including the effect of tails in the ra- 
diation reaction simply amounts to replacing the source 
Newtonian quadrupole moment lij by the tail-corrected 
expression [72] 
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where tq denotes some arbitrary time scale, for instance 



t-t' 
2Tn 



the one which appears in Eqs. (2.4). Note the factor 
AGM/c^ in front of the tail integral which is twice the 
factor 2GM/c^ in front of the tail integrals in (2.4). This 
factor ensures the consistency between the work done by 
the radiation reaction force in the local source and the 
total energy flux radiated at infinity from the source [15]. 

Thus the radiation reaction force including the 4PN 
tails takes (in harmonic coordinates) the same form as 
in Eq. (B3) but with replaced by /'J-''. This shows 
that the previous Newtonian argument still holds for the 
2.5PN -I- 4PN radiation reaction force and that the effect 
on the precession equation is still in the form of some 
irrelevant total time derivative: 
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Hence our conclusion that the constant-magnitude spins 
are secularly constant up to 4PN order corresponding to 
1.5PN radiation-reaction order. 



Appendix C: Gravitational-wave polarizations 

We derive in this Appendix the two gravitational- wave 
polarizations. They are computed from the projection 
formulas (2.2), using the expression (5.6) [together with 
Eq. (5.7)] for Shjj^ . We adopt the convention shown 
in Fig. 2 for the polarization vectors. To shorten the 
result, we denote the projections of the polarisation basis 
{N,P,Q} onto the moving triad {n, A,£} by e.g. P„, 
Pa, Pi- With this notation, we have 
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The authors can provide on demand a file containing the 
results in Mathematica® input format. 
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